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Mass effects in muon and semileptonic b — > c decays 

Alexey Pak and Andrzej Czarnecki 
Department of Physics, University of Alberta 
Edmonton, AB T6G 2G7, Canada 

Quantum chromodynamics (QCD) effects in the semileptonic decay b — > civ are evaluated to 
the second order in the coupling constant, 0(a^), and to several orders in the expansion in quark 
masses, m c /mb. Corrections are calculated for the total decay rate as well as for the first two 
moments of the lepton energy and the hadron system energy distributions. Applied to the muon 
decay, they decrease its predicted rate by —0.43 ppm. 



PACS numbers: 13.35.Bv,13.25.Hw,12.38.Bx 
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Decays of heavy fermions are an abundant source of in- 
formation about fundamental interactions. Particularly 
important among them is the muon (p) decay. Insensitive 
to strong interactions, it can be very precisely described 
by the electroweak model. The experiment MuLan at the 
Paul Scherrer Institute will likely measure the rate of the 
muon decay with an uncertainty better than 1 ppm and 
thus improve the determination of the Fermi constant Gf 
that describes the strength of the charged-current weak 
interaction Along with the fine structure constant a 
and the Z-boson mass, Gf is one of the three pillars of 
electroweak Standard Model tests Q. 

In a separate effort, the TWIST experiment at TRI- 
UMF measures energy and angular distributions of 
positrons in the p + decay, testing the Standard Model 
and searching for new interactions, notably new bosons 
predicted by left-right symmetric models [U 0, [1[ • 

To match this experimentalprogress, both the rate [|| 
and the energy distribution ij| have been calculated in 
quantum electrodynamics (QED) with 0(a 2 ) accurac 



Two-loop weak corrections have also been calculated 
In the decay rate studies, the electron mass m e was as- 
sumed negligible in the already small 0(a 2 ) effects. 

Here we show that the finite m e effect decreases the 
muon decay rate by about half ppm, exceeding previ- 
ous estimates [lfj and approaching the expected MuLan 
precision. 

The final-state fermion mass effects are much larger 
in the heavy-quark decay b — > civ. Studied in B- 
factories and the Tevatron, this process provides informa- 
tion about the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix element V c b, as well as about parameters governing 
heavy-quark dynamics (see Q for an up to date review 
and references). Also in this case, theoretical studies at 
0(a 2 ) are complete only for a massless final-state quark 
11]. For the actual massive c-quark, Oja 2 ) ef fects are 
known in some special cases of kinematics [12l . Ii~3l . flij . 
So-called Brodsky-Lepage-Mackenzie (BLM) corrections 
[HI have been obtained for the width [l6[ and moments 
of the energy spectrum [T3, EH- Also some logarithms 
of the mass m c have been determined to all orders in a s 
[lj|. Most recently, Melnikov calculated numerically the 
m c effects for the width and the first two moments of the 
energy distribution of hadrons and of the charged lepton 
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FIG. 1: Example of an expansion of a double-scale integral. 
Thick, thin, and dashed lines correspond to rat, m c , and 
massless propagators. Loop momenta can be all hard (Taylor 
expansion in m c of (a)), and one or more soft ((b), (c), and 
(d)). Double line in (d) denotes a static propagator. 



produced in this decay [20( . In this paper we present cor- 
responding analytical results obtained as an expansion in 
powers and logarithms of p = m c /mi,. 

The construction of such mass expansion is illustrated 
with an 0(a s ) example in Fig. [TJ Real and virtual cor- 
rections are calculated together as cuts of the diagram 
[TJa). Depending on virtualities of momenta flowing into 
and through the charm quark lines, m c can be treated 
as small compared to those momenta, or else those mo- 
menta can be treated as small compared to the 6-quark 
mass in other lines of the diagram. The most interesting 
case is shown in Fig. [IJd), where two loop momenta are 
of order m c . This configuration generates odd powers of 
p and will be discussed in some detail below. 

At 0(otg) the number of diagrams is larger and the 
analysis of momentum scales more challenging, with as 
many as 11 regions in some diagrams, but it follows the 
general pattern outlined above. As a result, even the 
four-loop diagrams shown in Fig. [2] are calculated as a 
series in p with exact coefficients. 

The most challenging contributions arise when all loop 
momenta are hard (~ mi,). All charm-quark lines are 
Taylor-expanded in m c , generating a large number of in- 
tegrals with various powers of denominator factors. The 
method of Ref. [2l| is employed to reduce all integrals to 
a set of master integrals. Most are known [lcj . sufficing 
to reproduce the m c = limit as well as the 0(p 2 ) terms. 



2 




pansion in p starting with 

2S 7T 2 

X l = — - Y -(34 + 241np)p 2 + 16^V 



273 
2~ 



- ( =^ -361np + 721n 2 p + 87r 2 ) p 4 



(3) 



FIG. 2: Examples of diagrams contributing to the 0{a 2 s ) 
rections to the 6-quark decay rate. 



The second order correction X2 is a sum of finite, gauge 
invariant parts proportional to various color factors, 

X 2 = C F X A + C A X N + T R {X C + X H + N L X L ). (4) 



Unfortunately, the following term 0(p 4 ) requires further 
e-expansion of the most difficult integrals (calculations 
are conducted in D = 4 — 2e dimensions). 

An approach to evaluating the additional terms is illus- 
trated with the the first diagram in Fig. [21 normalized at 
p 2 = — 1. After integrating the lepton loop momentum, 

I(ai,a 2 , 0,3, 0,4, a 5 , ae, a 7 , as, a 9 ; x) 
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with Z?i = k\, D 2 = (fci 4 
D A = (h + k 2 + k 3 + p) 2 , 



p) 2 , D 3 = {h + k 2 + 
D 5 = O2 + h + V? 



F x, 

Dg, = k 2 , and Dg = 2k 2 p. 



Dq = (k 3 + pY + x, D 7 = k 2 , us 
The artificial scale x, introduced in some propagators, 
must be set equal one for the needed on-shell case. 

This double-scale topology has 40 master integrals, 
instead of 21 in the single-scale case. 40 differen- 
tial equations are set up by relating master integrals 
to their derivatives, using — d x I({oi}; x) = (055+ + 
ag6 + )/({di}; x) (operator i + increments exponent Oi). 
Fortunately, this system splits into independent subsys- 
tems of at most four equations. Results are expanded to 
the needed order in e [221 ]. and integration constants are 
hxed using the behavior near x — > 00 [ljj. Finally, the 
limit x — > 1 is taken. 

In addition to this all-hard case, there are regions with 
some hard and some soft momenta, factorizing into dia- 
grams with less than four loops. Among them the most 
difficult ones have three soft loop momenta and are eval- 
uated using Ref . @|1|2|. 

Up to the two-gluon order, QCD corrections to the 
process b — > civ are parameterized by 



r(6 civ) = T C F 



C F 



-X, 



x 2 + 



where Tq 



Gj„\V cb \*ml 
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is the tree-level massless result, 



and a s is normalized in the MS scheme at to;,. The tree- 
level mass-dependent decay rate is 



X = l-8p 2 -24p 4 lnp- 



(2) 



The one-gluon correction, known exactly [26], has an ex- 



Xl, Xc, and Xh denote contributions of o, b-, and Nl = 
3 species of massless quarks. SU(3) color factors are Tr = 

1 4 /~i q 

7j, l~-F — 31 — 

Using techniques described above, expansions of X A , 
Xn, Xc, Xl, and Xh are obtained through 0(p 7 ), suf- 
ficient for sub-percent accuracy of X 2 for the physical 
p w 0.3. These expressions being rather lengthy, Table [J 
lists numerical values of the coefficients. Xm and Xc are 
shown explicitly for the purpose of subsequent discussion, 
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These correction significantly depend on p near its real- 
istic values: A 2 (0.25) = -6.59, while A 2 (0.3) = -4.89. 

In addition to the decay rate, corrections to the first 
two moments in lepton energy Ei — Ei/rrib, and the 
hadronic system energy Eh — E^/mb, have been com- 
puted. The average is taken over the whole phase space 
of decay products. These moments are parameterized by 
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and similarly for the moments of Eh, described by coef- 
ficients Hj. Table U shows the second-order corrections 
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Finally, Table [T] lists Uc, the c-quark contribution to 
r(6 — > uiD), defined by analogy with X c of Q. This 
result is useful e.g. in b — > S7 studies [27| , 

Our results can be tested by comparing logarithms of 
the mass ratio with the values predicted by the renor- 
malization group analysis [H^ |. That study summed up 
some of the logarithmic effects to all orders in the cou- 
pling constant. When those results are expanded in a s , 
three terms can be tested in order a 2 : p 2 In 2 p, p 3 \np, 
and p 4 In 3 p. They are related to the presence in the 
tree-level decay width ((2|) of terms p 2 and p 4 lnp, and 
the absence of p 3 there. 

Ref. [l9[ traced the origin of those terms to operators 
that can be constructed from the b and c quark fields 
and determined anomalous dimensions of those opera- 
tors. The logs in terms p 2 were found to arise from the 
running of the c-quark mass, and those in p 4 - from a 
complicated mixing of a variety of dimension seven oper- 
ators. Our results fully agree with that analysis. 

However, terms p 3 In p turn out to work differently In 
[l9| they were attributed to the running of m c , with the 
resulting coefficient that disagrees with our result. We 
find that those terms originate from the four-quark op- 
erator h^T^c cT^hb. Here hb denotes the static field 
describing a slow quark b. As discussed in Ref. [l9| . 
this operator gives rise to terms p 3 in the tree-level de- 
cay in the case of a vector coupling (T^ = 7^), but 
does not contribute at the tree-level in the chiral case 
(r M = 7^(1 — 75)). However, we find that at 0(a s ) it has 
a finite matrix element as shown in Fig. []Jd) , responsible 
for the cubic mass term in the decay width 

At the next loop level, the effects of the coupling con- 
stant running and of the anomalous dimension of that op- 
erator generate terms a 2 p 3 In p. Including charm, Nt + 1 
quarks contribute to the coupling running. Denoting 
p = %C A - ±T R (N L + 1), one finds 

- 8(3 + y T * - U(J a = ^T R N L + b ^T R - ^C A , (6) 

in agreement with the coefficients of 7r 2 p 3 Inp in Eq. ([5]). 

The linear term — 5ir 2 p/4 in Eq. ([5]) is noteworthy. As 
will now be explained, it arises because the on-shell (pole) 
definition of the 6-quark mass has been used. Although 
not suitable for phenomenology, it has been adopted for 
the ease of comparisons with Ref. [2(| and simplicity of 
presentation. The linear term m q /mb arises from a q- 
quark (m q <C mb) loop inserted in the gluon propagator 
in Fig. [ljd). Thus it does not depend on the final-state 
quark mass and equally affects decays b — > civ and b — > 
u£v. Arising from the gluons with momentum 0{m q ) this 
effect becomes a problem for the perturbative analysis 
when q is light, m q < Aqcd- This illustrates how a 
linear Aqcd /nib correction appears even in the total b- 
quark width when the pole mass is used. If a short- 
distance mass definition is used, such as the MS mass 



[28| , such terms are absorbed into the lowest-order decay 
width and are absent in higher-orders of the a s expansion 
[29I . f30l . HH . Note the factor five in the coefficient of the 
linear term, related to the fifth power of nib in the width 
formula. 

Where the linear correction really counts is the muon 
decay, whose width is traditionally expressed using the 
muon pole mass m^. The previous study of the muon 
decay [l0| neglected the electron mass m e . Its effect 
on the decay rate was assumed to arise from terms 
(«) 2 ( m*) in 2 (mA ~ 1 5 . 10 -8 The theoretical er _ 

Kir J J \«V J 

ror was estimated by taking the coefficient of this term to 
be 24. Due to the overlooked linear correction the elec- 
tron mass effect turns out to be even larger, affecting the 
determination of the Fermi constant from the anticipated 
new measurement. 

The Fermi constant is determined from the measured 
muon lifetime using the relation [1 01 ] 
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where Aq describes effects of the finite electron mass 
and radiative corrections in the limit of the four-fermion 
contact interaction (m M M\y)- The latter have been 
known exactly in the first order in a, and in the limit of 
zero m e in order a 2 . Present result extends the knowl- 
edge of the a 2 correction to the case of the finite m e . The 
extra shift is Aq(m e ) ~ —0.43 • 10~ 6 , comparable with 
the expected precision of the MuLan result. 

Even though the coefficient of the 

(B.) 2 (m^\ ln 2 fi^ N ) term is only -11. less than 

half the value taken in |K$ for the purpose of an error 
estimate, the total electron mass effect, dominated by 
the linear term m e /m^, is larger than expected. 

For the decay b — > civ, we confirm numerical results 
of [20] . The flexible numerical method developed in that 
reference enables one to impose phase space cuts. Our 
approach is complementary. It provides analytical re- 
sults and gives insight into the small-mass region, im- 
portant for the muon decay. It also facilitates changes 
of the scheme used for the heavy quark masses. Finally, 
it reveals the structure of logarithms and highlights the 
relative importance of various operators, improving on 
the analysis of [l9j]. But is it possible to combine the 
analytical approach with cuts on the lepton energy? 

According to (2Cj , effects of cuts for the lepton energy 
moments can be modeled with the tree-level distribu- 
tion [32]. For the cut E\ > E cut = 1 GeV one finds 
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For the rate 



L^cut the error is only —4%, and less than —2% for L^cut) 
indicating smallness of the hard-gluon radiation. 

The relative impact of cuts on the hadron energy mo- 
ments (EJ^) depends on the order n of the moment 

only very weakly To within 1%, H^JH^ and 



^2 cut/ H-f' both equal L^^JL^/' (and significantly ex- 
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TABLE I: Coefficients of p 1 \d? p in expansions of the components of X2, Eq. (|4]), and of the moments (Ef' h ). The last line 
shows the analogue of Xc for the decay b — > ulv. 



ceed H^/H^). 

ilarly: H^/H™ 
0.1%. Thus the effect of cuts E cut < 1 GeV on QCD 
corrections to (E%) can be approximated by that on the 
rate. 

Our result for the correction to the width exceeds the 
earlier estimate based on an interpolation of special kine- 
matic results 1J]. The full mass dependence being now 



First-order corrections behave sim- 
t/^ 2) « 4°LM 0) within 



-"l.cut/ 



known, it is clear that the correction varies strongly as 
a function of the mass ratio and is close to vanishing 
near the physical value. This indicates cancellations that 
were not taken into account in [141 ]. In particular, the 
non-BLM correction turns out to be about 1.7(ev s /7r) 2 



[20|. While the full phenomenological analysis requires 
a fit of the rate and moments, we expect this correction 
to decrease the value of \V c b\ determined from inclusive 
b decays by about one percent, bringing it slightly closer 
to the exclusive-decay result. 
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